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Abstract. —The southern Chinese Osmunda Xmildei has been suggested to be an intersubgeneric 
hybrid, i.e., O. japonica (subgenus Osmunda) X O. angustifolia (subgenus Plenasium ) or O. 
japonica X O. vachellii (subgenus Plenasium). These interpretations were based on morphological, 
cytological, and/or chloroplast DNA data, yet the parents of the hybrid remained unclear. 
Molecular phylogenetic relationships inferred here from chloroplast rhcL sequences and three 

nuclear DNA markers show that O. Xmildei is most likely a hybrid between the paternal O. 
japonica and the maternal O. vachellii. 

Key Words. —EST, intersubgeneric hybrid, Osmunda japonica , Osmunda Xmildei , Osmunda 
vachellii, rbcL 


The genus Osmunda of the leptosporangiate fern family Osmundaceae has 
natural hybrids (Kato, 2009; Tsutsumi et al, 2011). One such is Osmunda 
Xmildei C.Chr. ( = 0. bipinnata Hook., a later homonym of O. bipinnata L.), 
which nearly became extinct in its known range in Hong Kong. However, this 
hybrid was recently found in Shenzhen, Guangdong, and Mt. Qiyun, Jiangxi 
(Zhang et al., 2008), and less than 10 individuals are known. It was also found 
in Zhangjiajie, Hunan (Y.-H. Yan, pers. comm.). It can propagate via spores in 
experimental conditions (J.-F. Yang, unpubl. data), but it is uncertain if the 
individuals were derived from spore propagations or from independent 
formations of the hybrid. Osmunda Xmildei is characterized by subcoria- 
ceous, bipinnate-bipinnatifid leaves with round, entire pinnules, and fertile 
pinnae inserted below the middle of the leaf. For the origin of O. Xmildei, 
two possibilities were proposed (Fig. 1). Based on karyological and 
morphological analyses, He et al. (2006) suggested that O. Xmildei is a 
hybrid of O. japonica Thunb. (subgenus Osmunda) and O. angustifolia Ching 
(subgenus Plenasium ). Zhang et al. (2008) observed the absence of 
chromosome pairings at meiosis and resulting abortive spores in O. Xmildei, 
and suggested that it is a sterile Fl hybrid, but argued that the parents were O. 


* Corresponding author: (e-mail: tsutsumi@kahaku.go.jp) 
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O. angustifolia (A) x O. japonica (J) —► O. * mildei (M) 

A 



Fig. 1 . Candidate parentage of Osniunda Xmildei. A: O. angustifolia : pinna (on left side) and part 
(on right side). J: O. japonica : pinnule and part. M: O. Xmildei : pinna and part. V: O, vachellii : 
pinna and part. Scale unit is cm. Two proposed crossings are illustrated at the top and bottom. 


japonica and O. vachellii Hook, (subgenus Plenasium ), because O. vachellii 
co-occurs with O. Xmildei , but O. angustifolia does not occur in some of the 
localities of O. Xmildei (Fig. 2). Gou et al. (2008) also proposed that O. 
vachellii is the maternal progenitor of O. Xmildei, based on inferences from 
chloroplast DNA sequence data. Under either parentage hypothesis O. 
Xmildei is likely an intersubgeneric hybrid between the subgenera Osmunda 
and Plenasium. 

There are three more known hybrids reported in Osmundaceae (Kato, 2009). 
Eastern North American O. Xruggii Tryon is O. regalis L. (subgenus Osmunda) 
X O. claytoniana L. (subgenus Claytosmunda ) (Tryon, 1940; Wagner et al., 
1978; Whetstone and Atkinson, 1993; Li and Haufler, 1994). Japanese O. 
Xnipponica Makino is O. japonica (subgenus Osmunda) X O. claytoniana 
(subgenus Claytosmunda) (Ito, 1964), but Sugimoto (1979) suggested that it is 
an intergeneric hybrid, i.e., O. japonica X Osmundastrum cinnamomeum (L.) 
C.Presl. Neither is given molecular evidence. Finally, Japanese O. xintermedia 
(Honda) Sugim. is O. japonica X O. lancea Thunb. (both in subgenus 
Osmunda ; Tagawa, 1959: Iwatsuki, 1995; Tatuno and Yoshida, 1966; Yatabe 
et al, 20091. 
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Fig. 2. Map showing distribution of Osmunda angustifolia (solid squares), 0. japonica (broken 
line), 0. Xmildei (asterisks) and 0. vachellii (solid line). Solid asterisk shows the place of 0. 
Xmildei sampled in this study. 


The evidence offered to support the parentage of O. Xmildei in previous 
studies is inconclusive, because it lacked nuclear sequence data. Combining 
nuclear and chloroplast DNA data are useful to identify a maternal and a paternal 
parent. This study examines the aforementioned alternative origins of the hybrid, 
based on chloroplast rbcL sequences and three nuclear DNA sequences. 


Materials and Methods 


Materials .—Samples of three species of subgenus Plenasium ( O . vachellii, O. 
banksiifolia (C.Presl) Kuhn, O. angustifolia), the putative intersubgeneric 
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hybrid O. Xmildei, one species of subgenus Claytosmunda ( O. claytoniana), 
Osmundastrum cinnamomeum, and Todea barbara (L.) T.Moore, along with 
three species of subgenus Osmunda (O. japonica, O. lancea, O. regalis), were 

collected in the field or botanical gardens. The sources of materials used are 
shown in Fable 1. 

Sequencing of chloroplast and nuclear DNA. —Leaf fragments were used for 
molecular analysis. DNA was extracted from fresh or silica-gel-dried material 
using a QIAGEN DNeasy Mini Kit (QIAGEN, Valencia, CA) following the 
manufacturer’s instruction. Three nuclear DNA markers (EST_L058, 
EST_L110, EST_L258) selected from the expressed sequence tag library 
developed by Yatabe et al. (2009) and the chloroplast locus rbcL were 
analyzed. Primers for amplification and sequencing, and detailed information 
on the three nuclear markers are shown in Table 2 and Table 3, respectively. 
PCR was performed using a Perkin-Elmer 9700 DNA thermal cycler (Applied 
Biosystems, Foster, GA) with Ex Taq DNA polymerase (TaKaRa Bio, Tokyo, 
Japan) and Ampdirect Plus (Shimadzu, Kyoto. Japan) in 35 denaturation. 
annealing, and elongation cycles (30 sec at 94 C; 30 sec at 50 C for the rbcL and 
58 C in the three nuclear markers; and 90 sec at 72 C) with a final elongation 
step (7 min at 72 C). The PCR products were purified with ExoSAP-IT (USB 
corporation, Cleveland, OH) following the manufacturer’s instructions. 
Sequencing was conducted using an ABI PRISM 3130x1 Genetic Analyzer 
(Applied Biosvstems). The raw sequence data were assembled using Seqman II 
(Dnastar, Madison, WI). For the sample of O. Xmildei, PCR of the nuclear DNA 
markers was performed with PrimeSTAR Max DNA polymerase (TaKaRa Bio) 
in 35 cycles (10 sec at 98 C, 5 sec at 55 C and 5 sec at 72 C). The PCR products 
were cloned using a pGEM-T Vector System I (Promega, Madison, WI) and at 
least 10 clones were sequenced. Minor variants from single clones, presumably 
sequencing errors, were observed. Therefore a consensus sequence was used 
for each allele type; the differences between each consensus sequence and the 
original clones are shown in Table 4. The assembled sequences were aligned 
by Clustal X program (Thompson et al., 1997) and then aligned manually. 

Molecular phylogeny. —Phylogenetic analyses were performed by maximum 
parsimony (MP) and Bayesian analysis. Registered sequences of Osmundaceae 
in Genbank were added into the analyses (see Table 1). Maximum parsimony 
(MP) inference was conducted with PAUP* 4.0bl0 (Swofford, 2002). The bases 
that could not been identified were treated as unknown (N), and gaps were 
treated as missing data. All characters were equally weighed and heuristic 
searches were conducted with 1000 random addition replicates involving TBR 
branch swapping. Bootstrap values were calculated from 1000 pseudorepli- 
cates, each with 100 random additions. For the Bayesian analvses, MrMod- 

md tr 

eltest 2.0 (Nylander, 2004) was used to determine the nucleotide substitution 
model. Bayesian searches were conducted by MCMC with two independent 
sets of four chains, each run for ten million generations, sampling every 100 
generations by MrBayes 3.1.2 (Huelsenbeck and Ronquist, 2001; Ronquist and 
Huelsenbeck, 2003). The nucleotide model selected was: rbcL, GTR + I + G; 
EST_L058, HKY + I; EST_L110. GTR + I; EST_L258, GTR + I. The program 
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Table 2. Primers used for the analyses of three nuclear markers. 


Marker 

Primer 

Primer sequence 

Reference 

EST L058 

EST_L058F 

EST_L058R 

ATA AGGTTTCGCCCTCG AAT 
TCTTGGAGTTGCGAGTTCAC 

Yatabe et al. 2009 
Yatabe et al. 2009 

EST.LllO 

EST.L110F (P108) 
ESI L110R (P107) 

CATTGCA CATCGGAGATGAT 
CACAGCTGAAACACCCTGAA 

Yatabe et al. 2009 
Yatabe et al. 2009 

EST_L258 

EST_L258F 

EST L258R 

TCATGGCGACTGTGAAGAAG 

CGCCCTTTGGATTTACGATA 

Yatabe et al. 2009 
Yatabe et al. 2009 


Tracer (Rambaut and Drummond, 2009) was used to check the runs had 
reached stationarity and effective sample size of all the parameters was high 
(>100). The first 2.5 million generations before sufficient stationary genera¬ 
tions were discarded as burn-in periods and the rest of trees were used to 
calculate posterior probabilities. Osmundastrum cinnamomeum and Todea 
barbara were used as outgroups (Yatabe et al., 1999; Metzgar et al., 2008). 


Results 

Bayesian and maximum parsimony analyses of each dataset produced 
congruent topologies (Bayesian consensus trees shown in Figs. 3-6). Maxi¬ 
mum parsimony analyses resulted in three shortest trees of a length of 146 
steps (Cl = 0.77, HI = 0.23, RI = 0.93) for chloroplast rbcL (1227 bp) with 101 
parsimony-informative characters, 78 shortest trees of a length of 47 steps (Cl 
= 0.89, HI = 0.11, RI = 0.88) for nuclear EST_L058 (198 bp) with 18 
parsimony-informative characters, ten shortest MP trees of a length of 125 
steps (Cl — 0.91, HI = 0.09, RI = 0.91) for nuclear EST_L110 (572 bp) with 57 
parsimony-informative characters, and two shortest trees of a length of 122 
steps (Cl = 0.89, HI — 0.12, RI = 0.94) for nuclear EST_L258 (361 bp) with 59 
parsimony-informative characters. 


Table 3, Total lengths of three nuclear markers, coding and non-coding regions in Osmunda 
japonica (J3 in Table 1), and identified genes with one of the highest E-value (< 0.001) obtained bv 
Blast search (blastn) using sequences of ESI libraries (Yatabe et al. 2009). 


Marker 

Accession 

no. 

Total length 
(bp) 

Coding re¬ 
gion (bp) 

Non-coding 
region (bp) 

Putative gene 
(species! 

GenBank hit 
accession no. 
(E-value) 

EST_L058 

FS993661 

198 

101 

97 

glycolate oxidase 
(gox) 

i; Zantedeschia 
aethiopica) 

AYl 73074 
(2 X 10“ 124 ) 

EST LI10 

FS993713 

548 

148 

400 

glvcerol-3- 
phosphate 
[Zea mays) 

EU964956 
(9 X 10 -77 ) 

EST_ L258 

FS993861 

317 

79 

238 

ribosomal 
protein LI 7 
{Castanea sativa) 

AF334838 
(9 X 10 -102 ) 
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Table 4. Allele types, numbers of clones used, and variations between a consensus sequence and 
the original clones for Osmunda Xmildei samples in three nuclear markers examined. 




Allele 

(Jen bank ac¬ 
cession No. 



Number of clones: 


Included in Identical to 
consensus consensus Differing bv 
seq. seq. one subst’n 

Differing by 

two subst'ns 

Differing by 
three subst'ns 

EST 

L058 

A 

AB672788 

11 

7 

4 

- 

— 



B 

AB672789 

4 

4 

* 

- 

- 

EST 

o 

hJ 

i 

A 

AB672826 

10 

4 

3 

1 

2 



B 

AB672827 

4 

2 

1 

• 

1 

EST, 

L258 

A 

AB672863 

6 

3 

3 





B 

AB672864 

7 

4 

1 

2 



The rbcL sequence of O. Xmildei is identical to that of O. vachellii and also 
was similar to those of O. banksiifolia and O. javanica Blume (Fig. 3). The 
Osmunda Xmildei sequence was more distantly related to O. angustifolia (but 
with low support), and very far from O. japonica and other species of subgenus 
Osmunda. 

In the three nuclear markers, the O. Xmildei sample has two distinct allele 
types (Figs. 4-6). One type (Type A) had the same sequence as some plants of 
O. japonica, while the other (Type B) had the same sequence as O. vachellii (in 
EST_L58 and L110 in Figs. 4 and 5) or a sequence very similar to it (in 
ESTJL258 in Fig. 6). In each of the three nuclear-marker trees, O. Xmildei was 
more closely related to O. vachellii than to O. angustifolia. 

Discussion 

Our trees constructed from the chloroplast gene and nuclear DNA sequences 
agree with previous trees in the monophyly of the three subgenera of Osmunda 
(Yatabe et al., 1999; Gou et al,, 2008; Metzgar et ah, 2008). The rbcL 
phylogenetic relationships of the subgenera are the same as Yatabe et al.’ s 
(1999) from the same gene and Metzgar et al.’ s (2008) from seven chloroplast 
loci including rbcL, and different from Gou et al. ’s (2008) rbcL relationships. 
The relationships deduced from the three nuclear EST markers are not 
consistent with each other, but the relationship of the EST_L110 agrees with 
the rbcL relationship in the subgenera Osmunda and Plenasium being sister to 
each other. 

All the phylogenetic trees inferred from the three nuclear EST markers show 
that Osmunda Xmildei has two distinct allele types, and one is identical to 
those of O. japonica, and the other formed a monophyletic clade with those 
of Plenasium species (Figs. 4-6). It is suggested that O. Xmildei is an 
intersubgeneric hybrid between the subgenera Osmunda and Plenasium. The 
Type B alleles of O. Xmildei have the same EST_L058, EST_L110 sequences 
and the closest EST_L258 sequences to those of O. vachellii, suggesting that O. 
Xmildei is most likely derived by hybridization of O. japonica and O. vachellii 
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* O. banksiifofia, Bgl 
- O.javanica 

O. vacheffii , Vgl 

O. xmildei M 
— O. angustifofia 


Plenasium 


1.00 


97 



O. cfaytoniana CLgl 

O. cfaytoniana CL2 
O. cfaytoniana CLg2 
O. cfaytoniana CL1 


Claytosmunda 


1 00 


100 


100 
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c 
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Todea barbara 

Osmundastrum cinnamomeum Cgl 


- Osmundastrum cinnamomeum 


- Osmundastrum cinnamomeum Cg2 


5 changes 


on chloroplast rbcL (1227 bp). Values above branches 
indicate posterior probabilities (>0.9) calculated by Bayesian analysis and those below branches 
indicate maximum parsimony bootstrap values (>60). Thick branches are highly supported 
(posterior probabilities p>0.95 and bootstrap values >90). Arrowhead indicates O. Xmildei. 
Abbreviations of materials follow Table 1. 


(Figs. 4—6). The chloroplast rbcL sequence of O. Xmildei is identical to that of 
O. vachellii (Fig. 3), suggesting that it is the maternal progenitor of O. Xmildei; 
hence O. japonica is paternal. This suggested parentage agrees with Zhang et 
al. (2008) and Gou et al. (2008), who suggested O. vachellii as the maternal 
parent, based on chloroplast DNA data and distributional data. Osmunda 
hanksiifolia is also very closely related to O. Xmildei, however, comparative 
morphology does not support that O. banksiifolia is a parent, because it has 
prominently dentate pinnae, whereas O. vachellii and O. Xmildei (and also O. 

japonica) are both distinct with entire or somewhat serrate pinnae or pinna- 
segments. 

This study analyzed a sample of O. Xmildei from Shenzhen, Guangdong 
(Fig. 2). It has very low spore viability and a very low offspring reproduction 
rate even in carefully controlled culture conditions (Zhang et al., 2008; J.-F. 
Yang et al., unpubl. data). Considering the low reproductive ability and a few 
isolated localities in southern China, it is possible that O. Xmildei is of 
multiple origins, although no molecular evidence is available. Osmunda 
Xruggii of eastern North America (Connecticut and Virginia, USA) is an 
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1 change 


Fig. 4. Bayesian consensus tree based on nuclear EST_L058 {198 bp). Figures above branches 
indicate posterior probabilities (>0,9) calculated by Bayesian analysis and those below branches 
indicate maximum parsimony bootstrap values (>60). Thick branches are highly supported 
(posterior probabilities >0.95 and bootstrap values >90). Abbreviations of materials follow 
Table 1. Arrowheads indicate allele types obtained from O. Xmildei. Numbers of clones of each 
allele type are in Table 4. 

intersubgeneric sterile hybrid derived from O. regalis (subgenus Osmunda) 
and O. claytoniana (subgenus Claytosmunda ; Tryon, 1940; Wagner et al., 1978; 
Whetstone and Atkinson, 1993; Li and Haufler, 1994). Like O. Xmildei, O. 
Xruggii grows together with the parents in a few localities (Wagner et al., 
1978). Wagner et al. (1978) described that transplants produced fertile pinnae, 
but produced aborted spores and only univalent chromosomes at meiosis, 
suggesting its high sterility. Osmunda xintermedia, which is widely 
distributed across Japan, is suggested to be an intrasubgeneric hybrid of 
multiple origins from O. japonica and O. lancea, and it is self-reproducible 
(Shimura, 1972; Yatabe et al., 2009). The differing levels of fertility and 
sterility in the three hybrids may reflect close or remote phylogenetic affinities 
[Yatabe et al., 1999; Metzgar et al., 2008). 

Osmunda japonica is distributed in eastern Asia extending west to the 
Himalayas and south to northern Vietnam (Kato, 2007). It is a likely parent for 
O. Xmildei, O. Xnipponica and O. xintermedia, all three of which occur 
within the distributional range of O. japonica. Osmunda vachellii, the other 
parent for O. Xmildei, also co-occurs with O. Xmildei (Fig. 2). From this 
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Fig. 5. Bayesian consensus tree based on nuclear EST_L110 (572 bp). Values above branches 
indicate posterior probabilities (>0.9) calculated by Bayesian analysis and those below branches 
indicate maximum parsimony bootstrap values (>60). Thick branches are highly supported 
(posterior probabilities >0.95 and bootstrap values >90). Abbreviations of materials follow 

Table 1. Arrowheads indicate allele types obtained from O. Xmildei. Numbers of clones of each 
allele type are in Table 4. 


distributional pattern, along with the distributions of O. Xruggii and its 
parents, we suggest that overlap of the parental species allowed the 
interspecific hybridization relatively recently. On the contrary, a fertile 

wf w fc/ 

tetraploid species of hybrid origin between O. japonica and O. regalis 
(subgenus Osmunda), occurs in northern Central Laos distant from the 
distribution ranges of both parents, and in northern Myanmar, distant from 
Central India where O. regalis occurs (Kato, 2007; Tsutsumi et al., 2011). 
Tsutsumi et al. (2011) suggested that the hybrid species arose when the 
distribution ranges of the parents overlapped, a pattern different from the 
current pattern. 
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